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ABSTRACT 

On the basis of the universal gas fraction in clusters of galaxies, we estimate that the effective 
thermal conductivity required to balance radiative cooling in the cores, where the gas temper- 
ature is 3-10 keV, is about one tenth of the Spitzer rate. This confirms that thermal conduction 
can be important for the energy balance provided that it is not highly suppressed by magnetic 
fields in the gas. We determine the global effective conductivity in a sample of 29 clusters 
using published X-ray data on the inferred cooling rates and show that most lie between one 
and one tenth of the Spitzer rate. More work on the profiles in cooling flow clusters is required 
to test the conduction hypothesis further. We examine the possibility that conduction operates 
during galaxy formation, and show that it provides a simple explanation for the upper-mass 
cutoff in galaxy masses. 
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1 INTRODUCTION 

Recent X-ray data on the cores of clusters of galaxies from XMM- 
Newton and Chandra indicate that some heat flux must be balanc- 
ing the radiative cooling losses of the gas (Peterson et al. 2001; 
Tamura et al. 2001; Fabian et al. 2001; Johnstone et al. 2002, 
Bohringer et al. 2002; McNamara et al. 2001; Nulsen et al. 2002). 
The source of heat is unclear, with the central radio source being 
one possibility (Churazov et al. 2002; Briiggen & Kaiser 2001; 
Reynolds et al. 2002) and conduction being another (Narayan & 
Medvedev 2001; Gruzinov 2002; Voigt et al. 2002). Here we ex- 
amine conduction and its consequences in more detail. 

The observations show that large temperature drops, by a fac- 
tor of three or more, are common close to the centre of clusters. 
Thermal conduction is therefore an obvious heat transport process. 
The question is one of the level of conductivity and whether it can 
be high enough to balance radiative losses, and also whether it can 
do so in a stable manner over the long time-scales required. For a 
highly ionized plasma such as the intracluster medium the maxi- 
mum rate is expected to be that calculated by Spitzer (1962), which 
we hereafter refer to as Spitzer conductivity with coefficient ks- 
Magnetic fields in the plasma should suppress the conductivity by 
some factor / so that k — fns- The value of / is not clear, but 
could be as high as one third for a tangled field (see Narayan & 
Medvedev 2001 and Malyshkin 2001 and references therein). Sup- 
pression should be very high across field lines and appears to be 
high across the very abrupt temperature drops known as cold fronts 
found in some clusters (Ettori & Fabian 2000; Markevich et al. 
2000; Vikhlinin et al. 2001). 

In Section 2 we show that the level of conductivity required to 
balance radiative cooling in the core of a cluster, where the radia- 
tive cooling time is less than the age of the Universe, is less than ks . 
This means that if conduction is unimpeded then it can balance ra- 



diation. We then compare this result with an effective conductivity 
from observations, determined by 

L(< r) 

KcS = 47rr 2 dT/dr' (1) 

where L is the total luminosity within radius r at which the tem- 
perature gradient dT/dr is inferred. The sample used is from the 
brightest 55 in the sky studied using ROSAT data by Peres et al. 
(1998), supplemented by some luminous clusters from the sample 
of Allen (2000). We choose objects with known central emission- 
line nebulosities, indicating probable large temperature drops in the 
gas. 

If thermal conduction (i.e. electron motion) is relatively unim- 
peded in a radial direction, this forces us to consider that ion motion 
might not be as restricted as normally assumed. In other words, sed- 
imentation of heavy ions is a possibility, with implications which 
we consider in Section 3. 

Finally we discuss the influence of conduction on galaxy for- 
mation, showing that it cannot balance cooling in small-to-medium 
mass galaxies but can have a significant effect on massive galaxies. 
It may therefore be responsible for the high-mass cutoff to galaxy 
masses. 



2 THE EFFECTIVE CONDUCTIVITY K eB 

Consider the gas in the core of a cluster. Within radius r gas is ra- 
diatively emitting energy at a rate L leading to a decreasing temper- 
ature unless heat is supplied at a similar rate, which for conductivity 
means L = 47rr 2 K^. Ignoring a small factor which depends on 
the density gradients within that radius, and approximating dT/dr 



2 



as aT/r (i.e. T cx r a ), the effective conductivity required is 

KeB 



STa 



(2) 



A is the radiative cooling function and n e is the r.m.s. electron den- 
sity within r. Now, the radiative cooling time is given by 



^cool 

where n is the gas density, and so 
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From the virial theorem, GM/r » kT/m, where m is the mean 
molecular weight, so for a gas mass fraction / gas we have 

/gasfc T 

Kcff = 2 • (5) 

am z Gat coo i 

As we shall see, for typical cluster parameters this is less than, but 
reasonably close to, ks ■ 

We next consider the effective conductivity required by data, 
k' bS , for a sample of cooling flow clusters drawn from the study of 
Peres et al. (1998). We use those objects for which optical emission 
lines have been detected from the central cluster galaxy, indicating 
that some gas has cooled out there. Peres et al. (1998) tabulate L(< 
r C ooi), where r coo i is the radius at which i coo i = ^H^ 1 = 13 Gyr, 
as well as the outer (wide-beam) temperature of each cluster. As- 
suming that dT/dr = OAT/r, which is appropriate for conduc- 
tion balancing bremsstrahlung cooling above ~2keV in a spherical, 
constant pressure cluster (Fabian et al. 1994), and is also a reason- 
able fit to observed temperature gradients, we plot in Fig. 1 
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The sample has been supplemented to a total of 29 by some massive 
clusters from the work of Allen (2000), where the integrated mass 
deposition rate within the cooling radius is again calculated from 
imaging data (we use IRAS 09104+4109, Abell 963, Zwicky 3146, 
Abell 1068, RXJ1347.5, Abell 1835, MS2137.3-2353, Abell 2390 
and MS 1455.0+223). All exhibit optical emission lines. Where the 
measured temperature of the cluster is greater than lOkeV, the value 
is taken from Allen et al. (2001). 

The bolometric luminosity emitted from within the cooling 
radius is calculated from 



L{< r coo i = M 

m 



(7) 



where a factor of 2 is used to take into account the change in grav- 
itational potential energy of the inflowing mass which is assumed 
in the estimates of Allen (2000). 

Fig. 1 shows that most cooling flow clusters have n' cS between 
ks and /ts/10. Conduction may therefore be important, as noted by 
Narayan & Medvedev (2001) and Gruzinov (2002). Our inclusion 
of the relation from Equation ^, however, shows that all clusters 
where the central cooling time is less than a Hubble time and which 
have the universal gas fraction must lie in this region too. It may, 
therefore, be a coincidence that k' cS for clusters is close to «s ■ 

Conversely, the results of Equation ^ indicate that clusters lie 
in an interesting regime, where conduction is potentially very ef- 
fective, if it operates. Simple global estimates like those plotted in 
Fig. 1, however, are poor diagnostics as to whether conduction is 
actually taking place. Detailed study of the temperature profiles is 
required to determine whether conduction is important or not. We 
have begun such tests using Chandra data on cluster cores (Voigt 




Temperature (10 K) 

Figure 1. re c ff determined from Equation [] using data from (stars) Peres 
et al. (1998) for clusters which have optical line emission noted for the 
central cluster galaxy, and from (solid dots) Allen (2000). Lines for Kg and 
one tenth res are shown, together with Equation pi (dashed), for the case 
/ g as = 0.17, a = 0.4. 



et al. 2002), which supports conduction being an important factor 
in the energy flow from about 100 kpc into the innermost 20 kpc. 
Within the smaller radius, cooling may dominate and lead to the 
observed optical nebulosities and star formation. 

Loeb (2002) has noted that if conduction operates at the 
Spitzer rate throughout clusters and beyond, they would leak their 
thermal energy into the surrounding intergalactic medium. There- 
fore Spitzer rate conduction may be restricted to a cluster core. 
Perhaps the formation of a cooling flow and consequent gas in- 
fall leads to magnetic field reconnection and reordering and to a 
more radial field (Soker & Sarazin 1990) which then allows con- 
duction to become important (a possibility mentioned by Bregman 
& David 1989). In other words the cooling flow produces the con- 
ditions which allow conduction to operate. 



3 IMPLICATIONS FOR SEDIMENTATION 

We note that if thermal conduction is uninhibited by magnetic fields 
then so may be sedimentation within the cluster core. This could 
occur if, for example, a significant fraction of the magnetic fields 
were radial. Basically this means that ions heavier than protons 
will accumulate towards the centre relative to hydrogen. Although 
iron may therefore accumulate (Fabian & Pringle 1977) the effect 
is more important for helium nuclei (Gilfanov & Sunyaev 1984). 
Calculations have recently been performed for a Navarro, Frenk & 
White (1997) potential by Qin & Wu (2000). The abundance gra- 
dients seen in iron, and some other elements, in some cluster cores 
(Fukazawa et al. 1994; Ezawa et al. 1997, De Grandi & Molendi 
2001; Sanders & Fabian 2002; Johnstone et al. 2002) may there- 
fore be partly accounted for. 

If helium does undergo sedimentation, then the inferred gas 
density in the central cooling regions is reduced. Any stars formed 
from cooled gas will be helium rich and so have short main- 
sequence lifetimes (Lynden-Bell 1967). Such a radical possibility 
would imply significant increases in the star formation rate asso- 
ciated with the observed distributed blue light (Johnstone, Fabian 



On conduction, cooling flows and galaxy formation 3 



Conduction dominates 



Radiation dominates 



50 

Radius (kpc) 



100 



Figure 2. Regime (above solid line) where conduction dominates and (be- 
low) where radiative cooling dominates. An isothermal potential has been 
assumed and lines of constant mass (from the virial theorem) are plotted 
(dotted lines). It explains why gas cannot cool sufficiently to make massive 
galaxies above a few times 10 12 Mq. Along the dashed line the radiative 
cooling time is 5 Gyr. 



& Nulsen 1987; Allen 1995; Cardiel et al. 1998; Crawford et al. 
1999). This will be studied further in later work. 

We note that if conduction is suppressed then so will be sed- 
imentation, but if conduction is balancing cooling then sedimenta- 
tion does not necessarily operate, unless the field lines are essen- 
tially radial. 



4 GALAXY FORMATION 

Most models for galaxy formation are based on gas falling into 
dark matter potential wells, being heated and then cooling radia- 
tively to form stars (e.g. White & Frenk 1991). For massive galax- 
ies the cooling process resembles a cooling flow since the gas does 
not cool immediately but forms a slowly cooling and settling at- 
mosphere. A problem arises then in that there is no simple upper 
mass limit to a galaxy, and a supermassive galaxy could in princi- 
ple form at the centre of a cluster as a result of a massive cooling 
flow. In practice, the models assume that visible stars do not form 
in potentials above some mass (e.g. Kauffmann et al. 1999). We 
now investigate the effect of conduction on this cooling process. 

Following the models of White & Frenk (1991) and Kauff- 
mann et al. (1999) we adopt a simple isothermal distribution for 
the dark matter and give the gas a similar density distribution; 



/gasMp 

87rromr 2 ' 



(8) 



where Mo is the total mass within radius ro. The luminosity within 
ro is then 



/gasM 2 A 

327rr-Qm 2 



(9) 



We have assumed that gas extends only in to a radius of two 
thirds ro in order that the luminosity is not dominated by an un- 
physical central cusp. (We justify this factor from the expression 
for the cooling rate in the models of Kauffmann et al. 1999 of 
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Figure 3. Total energy required to stem cooling flows. The equivalent ac- 
creted black hole masses, based on a radiative efficiency of 10 per cent, are 
indicated. Symbols as in Fig. 1. 



M = 4immr^dr c /dt, r c being the cooling radius at time t, which 
with L — | — kT gives the above value.) 



2 m 

Setting 



Mo 



fcTrp 
Gm 



(10) 



from the virial theorem, and balancing L by a conductive flux gives 
r2 .A / kTr \ 2 , „, 

(ID 
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A 2 dT 



dr 



Let dT/dr = aT/r, k = k T 5/2 and for bremsstrahlung A 
A T 1/2 then 



T 



This is approximately 



/ gas A fc 



8(47rr) 2 K m 4 G 2 a' 
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If cooling by metals is important (T < 10 7 K) then A = A m and 
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has been taken to 

1 _ — l-fcr-3.5 



V 30kpc 

These constraints are shown in Fig. 2. A L 
be 10 _23 ergcm 3 s _1 and ko = 5 x lCT 7 erg cm -1 s _1 K 
We parametrize the gas fraction here as / = 0.1; the fraction drops 
as stars form, or have formed at earlier stages. Above the solid line, 
conduction can overwhelm radiative cooling and below, cooling 
dominates. Dotted lines of constant mass from the virial theorem, 
as well as a representative galaxy formation line (dashed) where the 
radiative cooling time (from Equation ^ is 5 x 10 9 yr, are shown. 
Conduction can have a significant effect on the formation of galax- 
ies with a total mass exceeding a few times 10 12 Mq. It is also 
relevant to the work of Kauffmann et al. (1999) who assume that 
gas cooling in potential wells with a circular velocity exceeding 
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350 kms" 1 (which corresponds to a temperature of about 10 7 K) 
does not form visible stars. 

Conduction provides a simple explanation for the upper-mass 
cutoff in galaxy formation. Something must be suppressing contin- 
ued cooling in massive haloes and conduction close to the Spitzer 
rate should become effective in the required mass range. If the alter- 
native is heating by a central black hole, then it is not obvious why 
this becomes effective in massive galaxies, where the atmosphere 
is most extensive, and does not dominate in lower mass ones. 



5 DISCUSSION 

The level of conductivity required to balance radiative cooling at 
the cooling radius in most of the cores of 29 cooling flow clusters 
lies within an order of magnitude of the Spitzer rate. This need not 
indicate that conduction is operating but may merely be a coinci- 
dence, given the age of the Universe and the universal gas fraction. 
The detailed shape of the temperature profile and a comparison 
with conduction models may be required to test whether conduc- 
tion is a significant factor in suppressing cooling flows. Essentially, 
Fig. 1 just shows that it is energetically feasible for conduction to 
play a significant role in cluster cores. 

The total energy required to suppress a flow can be consider- 
able. We plot this energy, determined from L(< r coo i)t, in Fig. 3. 
In the diagram t is the age of the Universe, and may be reduced 
by a factor of a few, but the energy is nevertheless large. Whatever 
the heat source, it is to be among the major heat flows in the Uni- 
verse. If accretion onto the central black hole is responsible, with 
the energy flowing out in jets (see e.g. Churazov et al. 2002, and ref- 
erences therein), then assuming an efficiency of 0.1 of the accreted 
rest mass, the accumulated black hole masses lie between about 
10 9 and 10 10 Mq, varying inversely with the assumed efficiency. 
This is often taken as the radiative efficiency of accretion, but it is 
likely to be significantly smaller for the total jet power dissipated 
as heat into the surrounding gas. This may more than counterbal- 
ance the likely factor of a few overestimate for the cluster age (say 
5 Gyr rather than the 13 Gyr assumed). We therefore consider that 
the masses indicated are lower limits for this process. 

Determining the extent to which conduction operates in clus- 
ter cores from observation will be difficult. Cold fronts indicate that 
that it is suppressed at least locally in some situations. The reason- 
able agreement over a range of radii which we found in a study of 
the K e ff profile in two clusters (Voigt et al. 2002) is good supporting 
evidence for conduction operating in general. Work is in progress 
to study more clusters. Of great importance is the magnetic field 
structure and its connectedness. If the intracluster medium consists 
of many separate small magnetic structures then large scale con- 
duction would be inhibited except where reconnection allows struc- 
tures to merge. Faraday rotation studies indicate in some objects 
that the fields have coherence scales of a few kpc (Carilli & Taylor 
2002 and references therein). The temperature profiles of clusters 
so far appear to be smooth on the scales studied, although projec- 
tion effects are important, and observational exposures generally 
not long enough to probe subtle variations. In some objects (e.g. 
the Perseus cluster, Fabian et al. 2000; and the Centaurus cluster, 
Sanders & Fabian 2002) spiral-like X-ray structures are seen which 
could reflect changes in conductivity associated with magnetic field 
structure. 
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